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Abstract 

We present the first complete 1-loop diagrammatic calculation of the cross sec- 
tions for the neutral Higgs production processes e'^e~ — > Z^h^ and e~^e~ — > A^h^ 
in the minimal super symmetric standard model. We compare the results from the 
diagrammatic calculation with the corresponding ones of the simpler and compact 
effective potential approximation and discuss the typical size of the differences. 



*on leave of absence from Institute of Theoretical Physics, Warsaw University 







1 Introduction 



The Minimal Supersymmetric Standard Model (MSSM) predicts at least one light neutral 
scalar Higgs particle with mass below ~ 150 GeV. In order to experimentally detect 
possible signals of Higgs bosons and to trace back as far as possible the physical origin 
of a produced scalar particle, detailed studies for the decay and production processes of 
Higgs boson are required. 

As has been discovered several years ago [1-3], radiative corrections in the MSSM 
Higgs sector are large and have to be taken into account for phenomenological studies. 
Several methods have been developed to calculate the radiative corrections to the MSSM 
Higgs boson masses, production and decay rates. Three main approaches have been used 
for the calculation of 1-loop corrections: 

a) The Effective Potential Approach (EPA) 0: This method computes the dominant 
correction to the Higgs boson mass spectrum and coupling constants in a simple and 
fast way and is thus most suitable for numerical calculations, e.g. in Monte Carlo 
studies. In its practical realizations, it is based on the following approximations: 

- only the corrections to the 2-point Green's functions are included, and only 
the contributions from top/stop and bottom/sbottom are taken into account 

- the Green's functions are evaluated at momentum = 0, thus neglecting the 
momentum dependence. 

The EPA leads to expressions for the cross sections and decay branching ratios in 
an effective Born approximation, where the tree level quantities are replaced by the 
corrected Higgs boson masses and mixing angles. 

b) The Renormalization Group approach (RGE) 0]. This method also leads to an 
effective Born approximation in the formulae for cross sections. The 1-loop corrected 
Higgs boson masses and couplings are obtained by using RG equations. Large 
logarithmic terms can be resummed, but realistic mass spectra are not covered by 
this method since it relies on an effective SUSY scale. The momentum dependence 
of the self-energies and 3- and 4-point functions are neglected also in this approach. 

c) The diagrammatic calculation (Feynman Diagram Calculation, FDC) Q: The 
masses are calculated from the pole positions of the Higgs propagators, and the 
cross sections are obtained from the full set of 1-loop diagrams contributing to the 
amplitudes. The version of ref. @] is based on the on-shell renormalization scheme. 
It takes into account: 
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- the most general form of the MSSM lagrangian with soft breaking terms, 

- the virtual contributions from all the particles of the MSSM spectrum, 

- all 2-, 3- and 4-point Green's functions for a given process with Higgs particles, 

- the momentum dependence of the Green's functions, 

- the leading reducible diagrams of higher orders corrections. 

This method is technically most complicated, but also most accurate at the 1-loop 
level and can be used as the reference frame for simpler methods. Moreover, it 
allows to study the full parameter dependence of cross sections and decay rates. 

The experimental searches for Higgs bosons at LEPl and studies for the future 
searches at higher energies conventionally make use of the most compact effective 
potential approximation. The phenomenological study in ref. ^ for LEP2 and higher 
energies, was done in the FDC, but still not complete at 1-loop order. A check of the 
quality of the simpler approximations and estimates of their reliability in various energy 
and parameter ranges have also not yet been performed so far. 

In this paper we present the first complete 1-loop diagrammatic results for the cross 
sections for the neutral Higgs production processes e'^e^ Z^h^ and e^e~ A^h^. 
In addition, we compare the results from the complete diagrammatic calculation with 
the corresponding ones of the simpler and compact effective potential approximation and 
discuss the typical size of the differences. 

Recently some papers on the leading 2-loop corrections to the CP-even MSSM Higgs 
boson masses have been published 0. The main conclusion is that 2-loop corrections are 
also significant and tend to compensate partially the effects of 1-loop corrections. The 
calculations are based on the EPA and RG methods. Since one of the main emphasis 
of our study is to figure out the difference between complete and approximate results 
in a given order, we have not implemented the 2-loop terms. They would improve the 
1-loop FDC results in the same way as the approximations and thus do not influence the 
remaining differences which can only be obtained by an explicit diagrammatic calculation. 

In Section |^ we give a short description of the computational schemes for the masses 
and cross sections. This is followed in Section ^ by a discussion of the numerical results 
and a comparison of the various approximations. 

2 Outline of the calculations 

The tree level potential for the neutral MSSM Higgs bosons can be written as: 

= mlHl + mlHl + e,,{ml^H{Hi + H.c.) + \{g^ + g'\Hl - Hlf + \g\H^H^f (1) 
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Diagonalization of the mass matrices for the CP-even and the CP-odd scalars, follow- 
ing from the potential (|l]), leads to three physical particles: two CP-even Higgs bosons 
H^,h^ and one CP-odd Higgs boson , and defines their tree-level masses tjih , rrih and 
rriA , with mn > rnh , and the mixing angles a, (3. 

The 1-loop radiative corrections significantly modify the tree level relations between 
the masses and mixing angles. The way of calculating the radiative corrections in the 
EPA and FDC methods is briefly described as follows: 

In the EPA, the tree level potential V^^'^ is improved by adding the 1-loop terms 0: 

y (1) (Q2) = y{0) ^q2^ ^ ^y(l) (g2^ ^2) 

where 

At""W^)^^ E SUA.' (log ^- I) (3) 

quarks \ / 

squarks 

V^^\Q'^) is the tree level potential evaluated with couplings renormalized at the scale 
Q^, and Str denotes the supertrace over the third generation of quark and squark fields 
contributing to the generalized squared mass matrix A4^. 

The full 1-loop potential V^^^ is rediagonalized yielding the 1-loop corrected physical 
masses Mh , Mh and the effective mixing angle cte// . Thereby Ma and tan /3 are usually 
considered as the free input parameters. Explicit formulae for the 1-loop corrected CP- 
even Higgs mass matrices can be found in first paper of ref . [0 . The physical Higgs boson 
masses and 1-loop mixing angles are used in the Born formulae for the production cross 
sections and decay rates. 

In many experimental neutral MSSM Higgs searches a simplified version of the EPA 
( "e-approximation" ) was used in the analysis of the experimental data. In this approxi- 
mation only the leading correction coming from the top and stop quark loops to the Higgs 
mass matrix is included. In addition, the effects of the sfermion mixing are neglected. In 
this case the corrected Higgs mass matrix has the simple form: 

_ sin2^ / cot PMI + tan PM^ -M| - Ml , 

JVlcp-e.en ^ ^ -Ml - Ml t an /?M| + cot /?M1 + 6 ' ^ ^ 

where e summarizes the leading 1-loop corrections: 



v^TT^sin^/? 



This yields: 



MJ, . M±M±i ± JM±M)i±i! _ MlMl oos'^ 2, + i^(Ml - Ml) ,6) 
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_ -(Ml + M|) tan /3 

tan«e//- Ml + Mlt^n^(5-{l + i^n^f3)Ml ^ ^ 

The 1-loop corrected Higgs boson masses , and the effective mixing angle thus 
depend on only one extra SUSY parameter, namely the product M^^ M^^ of the scalar top 
mass eigenvalues, in addition to the tree level parameters Ma and tan (3. 

In the FDC approach the 1-loop physical Higgs boson masses are obtained as the pole 
positions of the dressed scalar propagators. In order to obtain accurate values of the 
masses, also the leading reducible terms of higher order have to be taken into account. 
Mjj and M| are given by the solution of the equation of the form: 



Re 







For the calculations of the cross sections we need the full set of scalar 2-point functions, and 
also the 3- and 4-point Green's functions are taken into account. In Figure |I]the diagrams 
contributing to the e"'"e~ —>■ Z^h^,A^h^ process are collected. Specially important is the 
diagram f) in Figure |l|, illustrating the -H^ mixing on the external scalar line, as this 
mixing corresponds in a good approximation to the effect of introducing the effective 
mixing angle ag// into the Born formulae. The formulae for the cross sections obtained 
in the FDC now differ from the Born expressions, because not only the effective values 
of the masses are corrected but also new form factors and momentum dependent effects 
are considered. The analytic formulae for the cross sections for the processes e^e~ — > 
Z^h^,A^h^, including all but the box and 1-loop e'^e~h^ vertex contributions, can be 
found in the first paper of ref. [Q. Effects of box contributions are discussed in the 
third paper of ref. In the present paper we include also the 1-loop contributions to 
the e"'"e~/i° vertex in the cross section calculation. Note that the diagrams h), i), j) of 
Figure |I| are not part of the phenomeno logical study in |^ . 



3 Results on Higgs masses and production cross sec- 
tions 

In this section we present the results for Z^h^ and A^h^ production from the complete 
diagrammatic 1-loop calculation (FDC) and discuss the quality of simpler approximations: 

- the effective potential approach (EPA), as described above, 

- the simplified EPA with sfermion mixing and bottom/sbottom loop neglected ("e- 
approximation" ), based on eqs. (@-0). 
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In order to give an orientation about the accuracy of the approximations we show com- 
parisons of the various methods in the different regions of the parameter space. 

From the theoretical point of view, the most convenient parameters for the Higgs 
sector are the mass of the CP-odd Higgs boson and the ratio tan/3 = ^. From the 
experimental point of view it is more natural to use the mass Mh of the lighter CP-even 
Higgs instead of the formal quantity tan /3. This, however, is technically more complicated 
for the calculational scheme. In the first part of this discussion we follow the conventional 
way of presentation and choose tan/3 as a free input parameter. In the second part we 
replace tan (3 in favor of the physical neutral Higgs mass Mh . This has the advantage 
to directly get the theoretical results for the production cross sections for the case that a 
Higgs boson is found and its mass has been determined. 

As a first step, we put together the predictions of the various methods for the physical 
mass Mh for given M4 and tan /3. Figure ^ is based on the parameters listed in Table [^. 
H is the parameter describing the Higgs doublet mixing in the MSSM superpotential. M2 



Parameter 




Ma 






M2 




At = Ah 


Value (GeV) 


175 


70 


1000 


300 


1000 


500 


1000 



Table 1: Parameters used for the numerical analysis. 



denotes the SU(2) gaugino mass parameter. For the U(l) gaugino mass we use the value 
Ml = I tan^ 6*1^/1/2, suggested by GUT constraints, n, M2 determine the chargino and 
neutralino sectors (for the detailed expressions see for example [0). Mgq, Mgi, At and 
Ah are the parameters entering the sfermion mass matrices. For simplicity we assume a 
common value Mgg for all generations of squarks, and a common Mgi for sleptons. For 
the scalar top mass matrix one ha^: 



(i - I sin^ Ow)M^ cos 2/5 + + M^^ -mt{^i cot /3 + At) 

-mt (/i cot (3 + At) I sin^ OwM^; cos 2(3 + + M^^ 



In Figure ^ only the mixing parameters have been turned off, i.e. At,Ah,fi ~ in 
order to have a situation where the "e-approximation" should be most suitable. Indeed, 
the "e-approximation" and the full EPA are almost identical in this case. 

^Actually in the calculation we use 6x6 sfermion mass matrices, including the possibility of intergen- 
erational mixing . Effects of this mixing appear for example when the mass parameters for different 
sfermion generations are not equal to each other. 
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We define the relative differences for the masses with respect to the FDC as follows: 

= • (10) 

As shown in Figure ^, the values for obtained by FDC and EPA are rather close 
and differ by less than 4%. FDC gives lower values of Mh than the approximations. 
This is a general result, valid for all the parameter choices, and connected with the 
negative contribution to Mh from the gaugino/higgsino sector which is not included in 
the EPA calculations. Figure ^ also shows that the "e-approximation" for differs more 
significantly from the FDC result. The reason is just the neglected mixing term. 

The results for the heavier CP-even scalar mass Mh from FDC and EPA are again 
in good agreement. The "e-approximation" for this case yields masses typically 15%-20% 
different from the FDC values. 

Figures ^, ^ display comparisons of the cross sections for the processes e~^e^ 
Z^h^, A^h^ for two center-of-mass energies: 205 GeV (Figure ^ and 500 GeV (Figure 
The parameters are the same as the ones listed in Table |^. In analogy to (p^) we define: 

^FDC _ EPA,e 

dazh = -pDC azh = (r{e e ^ Z h) 



a 



Zh 

-.FDC _ EPA,e 

da^f^ = aAh = (r{ee ^ h A ). (11) 

^Ah 

As one can see from both figures, the EPA predictions follow in general more closely 
the FDC results than the "e-approximation" . The numerical differences can reach 10-20% 
at 205 GeV and 30% at 500 GeV. They become more important with increasing energies, 
exceeding 40% at 1 TeV. Note, however, that in the region of large cross sections the EPA 
accuracy is better (20% up to 500 GeV). 

The formally large relative deviations appearing in azh for large tan/3 are not im- 
portant from a practical point of view since the cross section is extremely small in this 
region. 

The less accurate number obtained in the "e-approximation" can again be addressed 
to the neglected mixing effects in the scalar quark sector. 

To give a more global impression of the typical size of the differences between the EPA 
and FDC results, we have chosen 1000 random points (for each y/s value in Table H) from 
the hypercube in the MSSM parameters space with the following bounds: 

0.5 < tan/3 < 50 5 GeV < Ma < 150 GeV 

-500 GeV < < 500 GeV 200 GeV < M2 < 1000 GeV 

200 GeV < Msg < 1000 GeV 100 GeV < M^i < 300 GeV 

-Msg < A = A < Msg 
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We calculated the quantities |(5Mf^^| , \SM§^^\ , |(5crff^| and \Safl^^\ and averaged 
them arithmetically over all generated points of the parameter space. The results are 
summarized in Table |^. It shows that the predictions of both methods deviate in particular 
for azh- 



Quantity 


=205 GeV 


=500 GeV 


v/s =1 TeV 




2.7% 


2.7% 


2.7% 




2.2% 


2.2% 


2.2% 




22% 


29% 


32% 




5.5% 


9.5% 


15% 



Table 2: Differences between the EPA and FDC predictions averaged over a random 
sample of parameters. 

We have analyzed also the dependence of the differences between the EPA and the 
FDC predictions on the SUSY parameters: sfermion and gaugino masses, fi parameter 
and sfermion mixing parameters. In most cases the variation of those parameters does 
not have a large effect on the size of the differences between the EPA and FDC. The only 
exception is the gaugino mass, which is completely absent in the EPA. For a relatively 
heavy {Ma > 50 GeV) the increase of M2 causes a slow change of 5M^^^ , Sia^h^ , 
Sa^h^ (see Figure ||). This effect can be quite dramatic for small Ma ~ 5 — 20 GeV and 
M2 around 1 TeV, but such low Ma values are however already excluded by the LEPl 
measurements. 

We now turn to the more physical parametrization of the cross sections in terms 
of the two Higgs boson masses Ma and M^ . This parameterization is more clumsy 
in the calculations, but it has the advantage of physically well defined input quantities 
avoiding possible confusions from different renormalization schemes. Technically, tan (3 
is calculated for given M^ , M4 and the other input parameters. Varying M^ {Ma and 
other input quantities fixed) we thus obtain tan/3 and azh, ctah as functions of Mh . The 
predictions for the cross sections are, however, in general not unique because there can 
be two solutions for tan (3, thus yielding two branches for fixed values of Ma ■ 

The solutions for tan f3 are shown in Figure |^ for the input parameters of Table 
again for the FDC and the simpler approximations. As can be seen from the figure, the 
differences between the tan P values obtained in the EPA and FDC can reach 20%. Also 
significant differences can occur for the cross sections, as displayed in Figure ^ where the 
predictions of EPA and FDC for the azh and aAh are plotted as functions of Mh . The 
typical size of differences between the methods is 10-20% for y/s =205 GeV. 
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4 Summary 



We presented the results from a complete 1-loop diagrammatic calculation in the MSSM 
for the cross sections for neutral Higgs production via e~^e~ —>■ Z^h^ and e~^e~ —>■ A^h^, 
including box diagrams and e~^e~ vertex contributions which were not part of previous 
studies for Higgs searches in the literature. Comparisons between the FDC predictions 
with the simpler EPA approximation have shown that the EPA has an accuracy of typi- 
cally 10-20% in the parameter regions where the cross sections are large. The differences 
become bigger with increasing energy. They may be specially important for future high 
energy e^e~ colliders. The use of the physical input variables Ma, Mh avoids ambiguities 
from the definition of tan (3 in higher order, but the observed differences remain of the 
same size. For a better accuracy, the full FDC would be required. It will also be necessary 
to incorporate the leading 2nd order terms since they are at least as big as the non-leading 
1-loop contributions. 

The library of FORTRAN codes for the calculation of the 1-loop radiative corrections in 
the on-shell renormalization scheme to the MSSM neutral Higgs production and decay 
rates can be found at the URL address 

http: / / itpaxpl .physik. uni-karlsruhe. de/^rosiek/neutraLhiggs.htnil 
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Figure 1: Classes of Feynman diagrams contributing to the e^e~ — > Z^hP process in the 
FDC approach. The diagrams contributing to the e^e~ —>■ A^h^ process can be obtained 
by changing into A° on the external line and skipping the diagrams i), j). 
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Figure 2: Comparison of the physical masses obtained in the "e-approximation'' 
EPA and FDC. Parameters as given in Table |I|. 
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Figure 3: Comparison of the physical /i° masses obtained in the "e-approximation" , EPA 
and FDC. Mass parameters as given in Table |1], but mixing parameters ^,At,Ai, ^ 0. 



11 



0.8 r 




0.2 
0.1 8 
0.1 6 
0.1 4 
0.12 

0.1 
0.08 
0.06 
0.04 
0.02 




1 1 1 1 1 1 1 


1 ' ' ' 

^ 




^Ah - 




EP.A 




f^Ah 




FDc : 




f^Ah 



tarifi 



0.6 r 



0.4 



-0,2 - 



-0.4 ■ 



-0.6 L 



tari/S 




tan/3 



Figure 4: Comparison of the cross sections a{e~^e Z^h^, h^A^) obtained in the " 
approximation" , EPA and FDC. Parameters as given in Table |l|, center-of-mass energ 
^/s = 205 GeV. 
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Figure 5: Comparison of the cross sections a{e~^e Z^h^, h^A^) obtained in the 
approximation" , EPA and FDC Parameters defined as in Table |l], center-of-mass energ 
^/s = 500 GeV. 
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Figure 6: Comparison of the cross sections a{e~^e~ — > Z^h^ , h^A^) as a function of M2 
in the "e-approximation" , EPA and FDC. Parameters as given in Table ||. tan (3 = 2 in 
the left plots and tan/3 = 50 in the right plots. Center-of-mass energy a/s = 500 GeV. 
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Figure 7: Comparison of the tan (3 dependence on Mh in the "e-approximation" , EPA 
and FDC Parameters as given in Table |l|. 5 tan /3 is the relative difference from the FDC 
values. 
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Figure 8: Comparison of the cross sections a{e~^e~ Z^hP ^ hPAP) as a function of Mh in 
the "e-approximation" , EPA and FDC Parameters as given in Table |I], center-of-mass 
energy y/s = 205 GeV. 
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